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Abstract 

We have calculated the large-g expansion for the energy cumulants 
at the phase transition point in the two-dimensional g-state Potts model 
to the 23rd order in using the finite lattice method. The obtained 

series allow us to give very precise estimates of the cumulants for q > 4 
on the first order transition point. The result confirm us the correctness 
of the conjecture by Bhattacharya el al. on the asymptotic behavior of 
the cumulants for q — > 4+. 
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1 Introduction 



It is a challenging problem to determine the order of the phase transition in the 
physical system where the order changes from the first order to the second one 
by a small variation of a parameter of the system. It is important to know how 
the quantities that remain finite at the first order transition point diverge or 
vanish when the parameter approaches the point at which the transition changes 
to the second order one. The g-state Potts model[[l], || in two dimensions gives 
a good place to investigate the properties of the phase transition around such 
a point. It exhibits the first order phase transition for q > 4 and the second 
order one for q < 4. Many quantities of this model at the phase transition 
point (3 = fit are known exactly for q > 4, including the latent heat|3| and 
the correlation length j|, ||, |(|, which vanishes and diverges, respectively, in the 
limit q — ► 4 + . One the other hand physically important quantities such as 
the specific heat and the magnetic susceptibility at the transition point, which 
also increase to infinity as q — ► 4 + , are not solved exactly. Bhattacharya et 
a/. |?J made a stimulating conjecture on the asymptotic behavior of the energy 
cumulants including the specific heat at the first order transition point; the 
relation between the energy cumulants and the correlation in their asymptotic 
behavior at the first order transition point for q — * 4+ will be the same as their 
relation in the second order phase transition for q = 4 and (3 — > fit, which is 
well known from their critical exponents. If this conjecture is true not only in 
the q-state Potts model but also in other general physical systems with a phase 
transition whose order change from the first order to the second one when some 
parameter of the system is varied, it would give a good guide in investigating 
such systems. 

Bhattacharya et aZ.|| also made the large-q expansion of the energy cumu- 
lants to order 10 in z = l/^/q and they analyzed the series according to the 
conjecture, giving the estimates of the cumulants at the first order transition 
point for q > 7. The obtained estimates are at least one order of magnitude bet- 
ter than those given by other methods such as the Monte Carlo simulations || 
and the low-(and high-)temperature expansions |TT| , in spite of the fact 
that the maximum size of the diagrams taken into account in their large-g ex- 
pansion is less than the correlation length for q < 15 while the Monte Carlo 
simulations are done carefully on the lattice with its size about 30 times as 
large as the correlation length and also the fact that the maximum size of the 
diagrams in their large-g expansion is as small as the half of the maximum size of 
the diagrams in the low-(and high-)temperature expansions These strongly sug- 
gest the correctness of the conjecture. The large-g series obtained by them are, 
however, not long enough to investigate the behavior of the energy cumulants 
for q closer to 4. 

In this paper we will enlarge the large-q series for the energy cumulants at the 
transition point to order of 21 or 23 in z using the finite lattice method |l3[ fbi] , 
R3, instead of the standard diagrammatic method used by Bhattacharya et al.. 



2 



The finite lattice method can in general give longer series than those generated 
by the diagrammatic method especially in lower spatial dimensions. In the 
diagrammatic method, one has to list up all the relevant diagrams and count the 
number they appear. In the finite lattice method we can skip this job and reduce 
the main work to the calculation of the expansion of the partition function for a 
series of finite size lattices, which can be done using the straightforward site-by- 
site integration jl6| |l7j without the diagrammatic technique. This method has 
been used mainly to generate the low- and high-temperature series in statistical 
systems and the strong coupling series in lattice gauge theory. One of the 
purposes of this paper is to demonstrate that this method is also applicable to 
the series expansion with respect to the parameter other than the temperature 
or the coupling constant. The long largc-g series obtained here by the finite 
lattice method enables us to examine the conjecture by Bhattacharya et al. on 
the asymptotic behavior of the energy cumulants for q very close to 4. We can 
also give the estimates of the energy cumulants for each q = 5, 6, 7, • • ■ that are 
orders of magnitude more precise than the estimates by the previous works and 
they would serve as a target in investigating this model in various contexts, for 
example in testing the efficiency of a new algorithm of the numerical simulation. 
A brief note on this subject was already presented in Jljfl , where we concentrated 
only on the specific heat. 

The remainder of this paper is organized as follows. The algorithm of the 
finite lattice method to generate the large-g expansion of the free energy are 
given for the disordered phase in section 2 and for the ordered phase in sec- 
tion 3. The relation between the two algorithm is also discussed in section 3. 
The obtained series for the energy cumulants at the phase transition points are 
presented in section 4. In section 5 the series are analyzed to investigate the 
asymptotic behavior of the cumulants at q — > 4 + as well as to make the estima- 
tion of the cumulants at various values of q > 4. In section 6 we conclude with 
a brief summary and a few remarks. 



2 Algorithm in the disordered phase 

The model is defined by the partition function 

Z=^exp // Y,*- - ■ W 

where represents the pair of nearest neighbor sites and Si takes the values 
1, 2, • • • , q. The phase transition point f3 t is given by exp (f3 t ) — I = ^fq. We first 
consider the free energy density in the disordered phase, which is given by 

F d {p)= lim {L x L y )- x \nZ d (L x ,L y ) . (2) 

Here the partition function Zd{L x , L y ) in the disordered phase should be cal- 
culated for the L x x L y lattice with the free boundary condition as depicted in 
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Fig.l. The large-g expansion of the partition function in this phase can be given 
through the Fortuin-Kasteleyn representation|l!| as 



with 



Zd{L x ., L y 



4>{X) 



-.L X L V 



XCA d 



(3) 



(eP - 1) 



y b(X) z 2[L m L y 



q c(X)-L x L y 



-c(X)]-b(X) 



(4) 



Here is the set consisting of all bonds in the lattice, X is any subset of 
Ad and represents a configuration of bonds, excluding the empty set (the null 
configuration) which gives the leading contribution and is separated in the first 
term of Eq.(§). In Eq.(@) Y = (e f3 - l)/y/q, z = 1/y/q, b(X) is the number of 
bonds in X and c(X) is the number of clusters of sites in X . Two sites connected 
to each other by a bond belong to the same cluster (An isolated single site is 
a cluster). The <f>(X) is called the activity of X. It is useful to note that the 
activity can be rewritten as 



cj>{X) 



O0 



l \K x \-b{x)+i(x) 



Y b(X) z b(X)-2l(X) 



(5) 



where l(X) is the number of independent closed loops formed by the bonds in 
X. 

We define Hd(l x , ly) for each l x x l y lattice (1 < l x < L x , 1 < l y < L y ) as 



H d (l x , ly) = In [Zd(l x , l y )/q 



(G) 



where Zd{l x ,ly) is the partition function for the l x x l y lattice, again with the 
free boundary condition as shown in Fig.l. We next define Wd{l x , l y ) recursively 
as 



W d {l x ,ly) =H d (l x ,l y ) 

- £ (l x -l' x + l)(l y -l' y + l)Wd(l' x ,l' y ) . (7) 

Here the term with l' x — l x and l' y — l y should be excluded in the summation. 
Then the free energy density defined by Eq.(Q) is given byjl5[ 

Fd(l3)=ln(q) + Y,W d (l x ,ly) ■ (8) 
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Now we have to know from what order in z the Wd(l x ,lp) starts. Only 
for this purpose, let us suppose to apply the standard cluster expansion^, pl| 
based on the diagrams the large-q expansion for the free energy in the disordered 
phase. It should be stressed that in the practical calculation of the finite lattice 
method we need not perform this cluster expansion at all. A polymer is defined 
as a set X of bonds which satisfies the condition that there is no division of X 
into two sets Xi and X 2 so that (f>(X) = 0(X!)</>(A 2 ). A set X composed of a 
single bond is a polymer. Below in this paragraph we discuss on the polymers 
that consists of two or more bonds. Now we define a cluster of bonds in X 
so that two bonds that are connected by a site belong to the same cluster. 
A polymer consisting of two or more bonds is composed of a single cluster of 
bonds, since a set of bonds composed of more than one disconnected clusters 
of bonds apparently does not satisfy the condition of a polymer. A cluster of 
bonds is a polymer if and only if the cluster of the bonds cannot be divided into 
tow or more clusters of bonds by cutting at any site that is shared by the bonds 
belonging to the cluster. We call such a cluster of bonds as one-site irreducible. 
In fact for a set X r consisting of a one-site reducible cluster of bonds that can 
be divided into clusters X\, X2, • ■ ■ , X n of bonds by cutting at a site, we have 
<p(X r ) = (j){Xi)4>(X2) ■ ■ ■ 4>(X n ), which implies that this set of bonds is not a 
polymer. For a cluster of bonds to be one-site irreducible, it is necessary that 
the bonds should form one or more closed loops and that each bond should 
belong to at least one of the loops, which we call closed. An example of a 
closed cluster of bonds is shown in Fig. 2(a). For each closed cluster of bonds, 
the loop of the bonds which forms the perimeter of the whole cluster is called 
the exterior boundary. The exterior boundary for the closed cluster of bonds in 
Fig. 2(a) is indicated in Fig. 2(b). An interior boundary is a loop inside which 
all bonds are removed and we call the part of the cluster surrounded by the 
interior boundary as a hole. The closed cluster of bonds in Fig. 2(a) has a hole 
with the interior boundary, as indicated in Fig. 2(c). For a cluster of bonds to be 
one-site irreducible, not only should it be closed but also its exterior boundary 
has no crossing. If the exterior boundary has a crossing, as in Fig. 3, the cluster 
of bonds is one-site reducible and is not a polymer. Thus a polymer is the closed 
cluster of bonds with its exterior boundary having no crossing. We list as an 
illustration in Fig. 4 all the polymers that can be embedded within 3x3 lattice 
with their activities. The polymers are skipped that can be obtained by rotating 
and/or reflecting those in Fig. 4. For later convenience we add some definitions 
concerning a polymer. A concave corner in the boundary of a polymer is the part 
of the boundary composed of two neighboring bonds that are bent perpendicular 
to each other toward the interior of the polymer. If two neighboring corners on 
the boundary are both concave, we call the two corners together with the bonds 
in a straight line connecting the two corners as a concave part of the boundary. 
Examples of the concave parts are given in Fig. 5. 

According to the theorem of the standard cluster expansion |20[ BjJ, the free 
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energy for the l x x l y lattice is given by 
\n.Z d {l x ,l y ) = lj y ]nq 

+ J2 c(X u X 2 ,---,X n )ct>(X 1 ) ( t>(X 2 )---4,(X n ), (9) 

{Xi .X2 ,•■ ■ ,X n -.connected} 

where the summation should be taken for every cluster of connected polymers 
{Xi, X2, • • • , X n } defined by the condition that there is no division of the cluster 
of polymers {Xi, X 2 , • • • , X n } into two clusters of polymers {X[,X 2l • • • , X' n } 
and {X[',X 2 ',- ■ -,X^_ m } (m<n) so that (j)(X 1 UX 2 U ■■■ UX n ) = <j>{X' 1 UX 2 U 
• • • U X' m )4>(X'{ UI^'U • • • U X%_ m ). (The 'cluster' of polymers here should be 
distinguished from the 'cluster' of sites in the Fortuin-Kasteleyn representation 
or the 'cluster' of bonds defined in the explanation of the polymers above.) The 
c(Xi, X 2 , • • ■ , X n ) in Eq.(||) is the numerical factor that depends on the way how 
the clusters of polymers are connected, the details of which are not necessary 
here. 

When a cluster of connected polymers can be embedded into the l x x l y 
lattice but cannot be embedded into any of its rectangular sub-lattices, we call 
it shortly to fit into the l x x l y lattice. Then we can prove the theorem that 
the Taylor expansion of the Wd(l x , l y ) with respect to z includes the contribution 
from all the clusters of polymers in the standard cluster expansion that fit into 
the l x x l y lattice. 

Now, according to this theorem, we can know the order of Wd(l x ,l y ) in z 
by examining the clusters of connected polymers that fit into the l x x l y lattice. 
We first note that Wd(l, 1) = since the lxl lattice consists of a single 
site with no bond, so that no polymer can be embedded into this lattice, and 
Wd(2, 1) = Wd(l,2) = 0{z) since the 2 x 1 or 1 x 2 lattice consists of two 
sites connected by one bond, so that a polymer composed of a single bond can 
be embedded into this lattice. We also note that Wd(l x , 1) = Wd{l, l y ) = for 
l x > 3 or l y > 3, because the corresponding lattice consists of the sites connected 
by bonds all on a straight line and the only clusters of polymers that fit into 
the lattice are those with each polymer composed only of a single bond and 
these polymers are not connected according to the definition of the connected 
polymers. 

We next consider the order of Wd(l x ,l y ) for l x > 2 and l y > 2. We first 
consider a cluster consisting of a single polymer that fits into the l x x l y lattice. If 
the polymer has a concave part in its exterior boundary we can add some bonds 
to make the part flat, obtaining a new polymer whose activity is changed by a 
factor of Y n z~ 1 with n the number of added bonds, so its order in z is lowered. 
This implies that a polymer that has concave part in the exterior boundary does 
not contribute to the lowest order term of Wd{l x , ly) in z. Let us next consider 
the case when the polymer has some holes. The interior boundary surrounding 
each hole has at least one concave part. (We take the hole as the exterior of 
the polymer when we say 'concave' here.) We can add some bonds to make the 
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concave part flat, again obtaining a new polymer with the order of its activity 
in z lowered by a factor of z^ 1 . This implies that a polymer that has a hole does 
not contribute to the lowest order term of Wd(l x ,ly) in z. Thus the remaining 
polymers that have neither a concave part on the exterior boundary nor a hole in 
the interior give the contribution to the lowest order term of Wd(l x , ly), examples 
of which are shown in Fig. 6. They can be converted into each other by adding or 
removing pairs of bonds to their concave corners or from their convex corners, 
respectively. When each pair of bonds is added or removed the activity of the 
polymer is changed by Y 2 with its order in z unchanged. Hence we see that 
their orders are all the same; that is, the order of z x+ v ~ 2 . Next we consider a 
cluster of connected tow or more polymers that fits into the lattice. It is enough 
to take into account the clusters of polymers that do not have either a concave 
part on the boundary or a hole in the interior, since we are going to know the 
clusters of connected polymers that contributes to the lowest order term in z. 
Among the clusters of connected two or three polymers that fit into the lattice, 
those give the lowest order contribution in which the polymers touch each other 
only through their exterior boundary and each neighboring pair of polymers 
share only one bond, like Fig. 7(a), and this type of clusters have the order of 
z ix+i y +n-3 f or cms ter of n polymers (n — 2,3), which is higher than the 
lowest order contribution from the clusters composed of a single polymer. We 
note that a cluster of two or three polymers in which each neighboring pair of 
polymers share only one site, an example of which is given in Fig. 7(b), is not 
connected according to the definition of the connected polymers. Among the 
cluster of connected more than four or more polymers that fit into the lattice, 
those give the lowest order contribution in which four of the polymers surround 
a plaquette and each neighboring pair of polymers share only one site, like Fig. 8, 
and this type of clusters have the order of z l *+ l y+ n ~ A ( n > 4) for the cluster 
of n polymers, which is also higher than the lowest order contribution from the 
clusters composed of a single polymer. In summary the lowest order term of 
Wd(l x ,ly) has the order of z l ^+ l v- 2 fo r l x > 2 and l y > 2, which comes from 
the clusters of a single polymer without a hole or a concave part fitting into 
l x x l y lattice. Thus in order to obtain the series to order z N for the free energy 
density, we have only to take into account in Eq.(||) all the rectangular lattices 
that satisfy l x + l y — 2 < N. 

We summarize the procedure of the finite lattice method to generate the 
large-g expansion series of the free energy density in the disordered phase. 

1) Calculate the partition function Zd(l x ,L y ) defined by Eq.(§) as a power 
series of z to order z N (and of Y in full order) for each l x x l y lattice with 
l x + l y — 2 < N. In these calculations use the transfer matrix method and inte- 
grate bond by bondjl^, 0, which reduces the necessary CPU time and storage. 

2) Calculate Hd{l x ,l y ) defined by Eq.(||) as a power series of z by Taylor ex- 
panding In Zd(l x , l y ) for each l x x l y lattice. 

3) Calculate Wd(l x ,l y ) defined recursively by Eq.(Q) for each l x x lattice. 

4) Sum up all the Wd(l x , lyVs with l x + l y — 2 < N according to Eq.rt|) to obtain 
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the expansion series for the free energy density. 

3 Algorithm in the ordered phase 

Next we consider the large-g expansion of the free energy density F ((3) in the 
ordered phase, which is given by the same equation as Eq.(|^) with the suffix 
d replaced by o where the partition function Z {L x ,L y ) in the ordered phase 
should be calculated for the L x x L y lattice with the fixed boundary condition 
as depicted in Fig. 9. The values of the spins on the boundary are all fixed to a 
single value among 1, 2, • • • , q. 

By the duality the ordered and disordered free energy densities are related 

as 

W)=F ( / 8) + 21n((e^-l)/Vg) 

for [f? - l)(e^ - 1) = q . (10) 

(n) 

One method to obtain the large-q series for the energy cumulant F at p = pt 
in the ordered phase is to use the relations 
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-(z - 90z 3 + 360z 5 )[i^ 1} - F G (1) ] , 

(11) 

which are derived directly from the duality relation (|Io|). These were given in 
Ref.]8| except for the last one. 

In spite of that, we give below the algorithm to generate the series for the free 
energy density in the ordered phase by the finite lattice method without relying 
on the duality relation. This algorithm will prove important when we are going 
to extend our work to generate the magnetization cumulants for instance, since 
they can be obtained by the derivative of the free energy density with respect 
to the external magnetic field and we have no longer duality relation for the 
free energy density in the presence of the external magnetic field. The large-g 
expansion of the partition function in this phase can be given also through the 
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Fortuin-Kasteleyn representation as 



Z (L x ,L y ) = (eP-l) 



X 



(12) 



with 



^(X) = (e^-l)- 6 ^^) 

= y KX) z b{X)-2c(X) _ 



(13) 



Here X is a configuration of the removed bonds and we also use X to represent 
the set of removed bonds. In the leading configuration for the fixed boundary 
condition all the pairs of nearest neighbor sites are connected by the bonds 
and there is no removed bonds (X is the empty set), and its contribution is 
separated as the first term of Eq.©. In Eq.(||) Y = *Jqj(e? - 1) and b(X) 
is the number of bonds in X and c(X) is the number of clusters of sites in the 
configuration of bonds with X removed without taking into account the clusters 
including the sites on the boundary. We define H Q (l x , l y ) for each l x x l y lattice 
(1 < l x < L x 



1 < l v < L y ) as 



H {l x ,l y )=\n Z (l x ,l y )/( e P -1) 



21. r 



(14) 



where Z (l x ,l y ) is the partition function for the l x x l y lattice, again with the 
fixed boundary condition as shown in Fig. 9. We next define W (l x ,l y ) recur- 
sively in the same way as in Eq. (R) . Then the free energy density in the ordered 
phase is given by 



F {p) = 2 In (e 13 - 1) + J2 W o{h, ly) 



(15) 



It is useful to compare now the large-g expansions by the finite lattice method 
for the disordered phase presented in the previous section and for the ordered 
phase presented here. We first note that the lattice with the free boundary 
condition in Fig.l and the lattice with the fixed boundary condition in Fig. 9 
are dual to each other in the sense that each site of one lattice is located at the 
center of a plaquette of the other lattice, when the two lattices have the same 
size. We denote X ~ X as the relation between the configuration X of the 
bonds on the lattice with the free boundary condition and the configuration X 
of the removed bonds on the lattice with the fixed boundary condition, if the 
location of each bond (identified with its center) in X coincides with the location 
of each corresponding removed bond (also identified with its center) in X (see 
Fig. 10(a) and (b)). Then, the number 1{X) of independent closed loops for each 
configuration X in Eq.(0) is equal to the number c(X) of clusters of sites in the 
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configuration with X removed in ( |l3| ) for X ~ X (see Fig. 10(a) and (b)), so 
that <f>(X) = ip(X) when Y = Y. We therefore know from Eq.(|), (§), @ and 
(P|) that Hd(l x ,l y ', Y, z) — H (l x , l y ; Y, z) when F = F, hence the same is true 
for WdQxi ly\ Y, z ) an d W (l X7 l y ;Y , z). Thus W D (l x , l y ) has the same order in z 
as Wd(Za;,Zy), that is, the order of z l * +ly ~ 2 . So in order to generate the series 
to order z N we have only to take into account the rectangular lattices with 
l x + l y — 2 < N in Eq.(p^). This conclusion can also be reached by repeating 
the discussion based on the standard cluster expansion using the diagrams. The 
discussion goes completely in parallel with those given for the disordered phase 
in the previous section if the set of bonds X representing the configuration of 
bonds for the disordered phase is replaced by the set of removed bonds X for 
the ordered phase. 

The relation of Wa(l x , l y \ Y, z) — W {l x ,l y ] Y, z) for Y = Y also implies that 
the free energy densities Fd((3) given by Eq.(||) and F ((3) given by Eq.(^5|) 
satisfy the duality relation (10), because the condition that Y = Y is equivalent 
to the condition for [3 and (3 in the duality relation in Eq.([To|). 



4 Series 

The algorithms in the previous two sections give the large-q expansion of the free 
energy density at the arbitrary inverse temperature (3. The energy cumulants 
can be obtained by taking its derivative with respect to (3. If we are interested 
only in the energy cumulants at the phase transition point, as is the case here, 
we can set Y = 1 + y and Y — 1 + y (with y — and y — at [3 — A)@ , then 
in order to obtain the n-th cumulant we have only to keep the expansion with 
respect to y to order y n as 



F. 



(n) 



df3~" 



F d {(3) 



P=p t 



E 



(16) 



using ^ = (1 + y + z)-£ and 



d n 

W 1 



F (f3) 



P=p t 



(17) 



using ^ = [l + y+il + yfz]-^. 

Using this procedure we have calculated the series to order N — 23 in z for 
the n-th order cumulants with n — 0, 1, 2 and to order N — 21 for the cumulants 
with n = 3, • • • , 6. The coefficients of the series are listed in Table 1.1 and 1.2 
for the disordered phase and in Table 2.1 and 2.2 for the ordered phase. We 
have applied the finite lattice method both in the ordered and disordered phases 
without relying on the duality relations ([n]) for the cumulants. We have checked 

+ l y — 2 < N have the correct 



that all of the Wd(l x ,l y ) and W (l x ,l y ) with l x 
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order in z as described in the previous two sections. The obtained series for the 
zcroth and first cumulants (i.e. the free energy and the internal energy) agree 
with the expansion of the exactly known expressions. The coefficients for 
agree with those by Bhattacharya et al. to order 10. We have also checked that 
the series satisfy the duality relations (|il"l ) for the cumulants. 



5 Analysis 

In this section we analyze the large-g series for the energy cumulants presented 
in the previous section. The asymptotic behavior of the energy cumulants at 
(3 — fit for q —>■ 4 + was conjectured by Bhattacharya et aZ.|?J as follows. It is 
well known that there is an asymptotic relation between the energy cumulants 
and the correlation length £ for (3 — > f3 t at q — 4 as 

(-1)«F W ~ A T ^~^ e n/2 - 2 , (18) 

with the constant A that is independent of n. This relation comes from the fact 
that the critical exponents of the correlation length and the specific heat in the 
second order phase transition are v — 2/3 and a = 2/3, respectively, for q = 4. 
Their conjecture claims that the relation ( |l8| ) is also held as the asymptotic 
relation between the energy cumulants and the correlation length at the first 
order phase transition point for q — > 4+. Then, from the known asymptotic 
behavior of the correlation length at the first order phase transition point for 
1 -» 4+ as|j 




where x = exp ( f^J with 8 defined by 2cosh# = ^fq (6 ~ y/q — 4/2 for q — > 4 + ), 

the conjecture implies that the asymptotic form of the energy cumulants in the 
limit of q — > 4 + is given by 

F W j (_l)»F o («) ^ ag n-2 r ^-|) a ,3n/2-2 _ 

We notice that the constants a and B should be independent of n and from the 
duality relations (|ll|) that they should be common to the ordered and disordered 
phases. 

To analyze the large-g series for these quantities, Bhattacharya et al. used a 
sophisticated method. Let f(z) be a function that has the asymptotic behavior 
of the right hand side of Eq. (Jl9|) . Then #ln/(z) is a regular function of z 
except for the possible terms proportional to the inverse powers of x, which 
vanishes rapidly when 9 approaches zero. Although 9 itself cannot be expanded 
as the power series of z since 9 diverge at z = 0, the quantity 1 — e~ e = 
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1 — 2z/(Vl — 4z 2 + 1) behaves like 8 for 9^0 and can be expanded as a 
power series of z. Then one may expect that (1 — e~ e ) ln/(z) will be a smooth 
function of z and it's Pade approximant will give a convergent result for q — > 4+. 
In fact using the 10th order series Bhattacharya et al. found that the Pade 
approximants of F^ 2 ' give convergent results for q > 7. The obtained results 
were more precise than the results of the Monte Carlo simulations || and the 
low- (or high-)temperature series jl(| O, [l2| as mentioned in the introduction 
and exhibits rather flat behavior of F^/x around 7 < q < 10, suggesting the 
correctness of the asymptotic behavior (|l^) for n = 2 with a ~ 0.76. If we 
use our longer 23rd order series of F^> to repeat these Pade analysis, we find 
that the convergent region of q is extended down to q ~ 5, with the accuracy 
of a few percent at q — 5 and the data of fx is rather flat for 5 < q < 10, 
which gives a stronger support to the conjecture that F^ 2 > is approaching it's 
asymptotic form jl9|). 

We, however, adopt here another method of the analysis that is simpler and 
more efficient. The latent heat C has the exact asymptotic form as|7| 

C ~ 3ttx- 1/2 , 

so, if fW has the asymptotic form of Eq.([n]), the product F^ C p is a smooth 
function of z for q — > 4+ when p = 3n — 4, while it behaves like some positive 
or negative powers of x and has an essential singularity with respect to z at 
q = 4 when p =/= 2>n— 4. Thus if we try the Pade approximation of F^C P it will 
lead to a convergent result for p — 3n — 4, while the approximants will scatter 
when p is off this value. In the Pade approximation of F^C the series for 
the ordered and disordered phases should be treated in a slightly different way, 
corresponding to the difference of the lowest order of the series; the series for 
the ordered phase starts from the order of z 2 , while the series for the disordered 
phase starts from the order of z. Thus we should solve the following equation: 

F („) £ 3n-4 = z rp M ( z yQ L ( z ) + 0(z M+L+r+1 ) , (20) 

where r is 2 for the ordered phase and 1 for the disordered phase, and Pm(z) 
and Ql{z) are polynomials of order M and L respectively with M + L + r = N 
where N is the truncated order of the series. Fig. 11 is the result of this analysis 

(2) 

of F^ C p for q = 4. We see very clearly that the Pade approximants have 
the best convergence at p = 2.00(1) and the convergence becomes bad rapidly 
when p leaves this value. All the other approximants for F^ n J C p have similar 
behaviors with respect to p with the best convergence at p = 3n — 4, although 
we do not show them in figures. These results are enough to convince us that 
F^ n J has the asymptotic behavior that is proportional to x 3n / 2 ~ 2 . 

Then we calculate the amplitude of F^ for each q by taking the Pade 
approximants of F ( n )£ 3 ™- 4 anc l multiplying £ 4_3n that is calculated from its 
exact expression. Usually a diagonal approximant (that is, M = L) or it's 
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neighborhood produces good estimates. We, however, adopt here all the ap- 
proximates with M > 8 and L > 8, since we have found in this case that each 
of these approximants is scattered almost randomly around their average value 
and we have found no reason to give a special position only to the diagonal 
approximant. Some denominators of the approximants have zero at z < 1/2 
(q > 4) and we have excluded such approximants. The results for n = 2, • • • , 6 
in the ordered and disordered phases are listed in Table 3-6 for integer values 
°f <7 > 5, with the second energy cumulant transformed to the specific heat by 
multiplying 01 . In Figs. 12-16, we also plot \F^ n '\/x~^'~ 2 versus 9 for n =2-6, 
respectively. In these figures the dotted and dashed lines indicate the errors for 
the disordered and ordered phases, respectively. In these tables and figures the 
given errors represent the magnitude of the fluctuation in the data of the Pade 
approximants. 

In any case, the estimates have surprisingly high accuracy. Our data are 
orders of magnitude more precise than the result of the Monte Carlo simulation 
by Janke and KapplerQ for q — 10, 15, 20 (also listed in Table 3-6), which was 
performed carefully in large enough lattices checking the finite size effects, and 
they are consistent with each other except for small differences in the data of 
the 5th and 6th cumulants for the ordered phase. On the other hand, the results 
of the low-temperature series [|l2| (also listed in Table 3 and 4) are inconsistent 
with these. This would be because we do not have the correct knowledge (or 
conjecture) of the possible singularity structure of the cumulants at some (3 < fit 
in extrapolating the low-temperature series toward (3 — (3 t - We should empha- 
size that the values of the specific heat has an accuracy of about 0.1 percent 
at q = 5 where the correlation length is as large as 2500 lattice spacing|(| and 
there is no previous data to be compared with, and of a few percent even at 
the limit of q = 4 where the correlation length diverges. The accuracy is better 
by a factor of more than 10 when compared with the result of using the Pade 
approximation for (1 — e~ e ) hiF^ n ', The reason for this is the following. Let us 
suppose that the Pade approximants gives estimates of F^C 3n ^ 4 — c(l ± S) 
and (1 - e- )\nF^ = c'(l ± 8'), then the estimate of from the former 
would also have an ambiguity factor of 1 ± S, while the estimate of F^ from 
the latter would have an ambiguity factor of e ±c s ^ e , which grows up for the 
plus sign and shrinks to zero for the minus sign when 9 — ► 0(q — > 4 + ). 

Next we determine the constant a and j3 in the asymptotic form (|w|) for 
q — > 4 + . The constants a can be determined by the value of F^ /x in the limit 
of q — » 4 + . We read from Fig. 12 that 



On the other hand, the constant B is given by the asymptotic value of the 
following combination for n > 3: 



a 



0.073 ± 0.002 . 




(21) 
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Their asymptotic values should be common to each n. The behavior of the 
quantity in Eq.(|2l|) for the ordered and disordered phases are given in Fig. 17 
and 18 respectively. The graphs for n = 3 and n = 4 show clear convergence 
even at 9 = giving the value of B very close to each other, that is, B ~ 0.38(3) 
We also notice that, although the results for n = 5 and n — 6 have considerably 
large errors in the small 6 region so we cannot say anything definite, the values 
of the quantity in Eq. ( ^l|) at the points where the errors begin to grow are also 
close to this value of 0.38(3). Thus we may estimate the constant B as 

B = 0.38 ±0.03 . 

In conclusion the Pade data of the large-g series convince us that the conjecture 
( |l9l) for the asymptotic form is true at least for n — 2, 3 and 4 and suggest its 
correctness for n — 5 and 6. 



6 Summary 

The finite lattice method were applied to generate the large-q expansion of the 
n-th energy cumulants [n = 0, 1, • • • , 6) at the phase transition point for the 
ordered and disordered phases in the two-dimensional g-state Potts model. It 
enabled us to extend the series to the 23rd order in z = 1/ yfq for n = 0, 1, 2 and 
to the 21st order for n = 3, • • • , 6 from the 10th order generated by Bhattacharya 
et al. using the standard diagrammatic method. We made the Pade analysis of 
the series for the product of each cumulant and some powers of the latent heat 
that is expected to be a smooth function of z if the conjecture by Bhattacharya 
et al. is true that the relation between the cumulants and the correlation length 
in the asymptotic behavior for /? — > (3t at q = 4 will be kept in their asymptotic 
behavior at the first order transition point for q — > 4 + . It allowed us to present 
very precise estimates of the cumulants at the transition point for q > 4. The 
accuracy of the estimates is orders of magnitude higher than the results of the 
Monte Carlo simulations and the low-(and high-) temperature series for q > 7, 
and for instance of about 0.1 percent at q — 5 where the correlation length 
is as large as a few thousand lattice spacing and there is no previous data to 
be compared with, and of a few percent even at the limit of q — 4 where the 
correlation length diverges. The resulting asymptotic behavior of the cumulants 
for q close to 4 is consistent with the conjecture not only for the second cumulant 
(or the specific heat) but also for the higher cumulants. 

Now it is quite natural to ask whether the conjecture by Bhattacharya et 
al. for the energy cumulants can be extended to other quantities such as the 
magnetization cumulants at the first order phase transition point that diverge 
in the limit of q — > 4 + . We can expect that the long series of the large-g 
expansion for the quantities would serve to answer this question as the long 
series for the energy cumulants did in this paper. It is rather straightforward 
to extend the large-q expansion by the finite lattice method to such quantities. 
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They are obtained in general by the derivative of the free energy density in the 
presence of appropriate external fields or sources, which can also be given by 
the algorithm of the finite lattice method described in this paper with a minor 
modification. Since no duality relation is known for the free energy density in 
the presence of the external fields, both of the algorithms presented here for 
the ordered and disordered phases will be useful, especially in what type of the 
boundary condition one should adopt for the finite size lattices. The calculation 
and analysis of the large-g series for the magnetization cumulants are now in 
progress. f23f 
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Figure 1: L x x L y lattice (in this figure L x = 4 and L y = 5) with the free 
boundary condition for the disordered phase. 
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(a) 



Figure 2: (a) Example of a closed cluster of bonds with a hole, (b) The same 
as (a) with the exterior boundary highlighted by thick lines, (c) The same as 
(a) with the interior boundary highlighted by thick lines. 



Figure 3: Example of a polymer whose exterior boundary has a crossing. 
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Figure 4: All the polymers that can be embedded within 3x3 lattice with 
their activities. 



Figure 5: Example of a polymer with two concave parts highlighted by thick 
lines. 
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Figure 6: Examples of the cluster consisting of a single polymer that contribute 
to the lowest order term of the Wa(l x , l y ) with l x = 4 and l y = 5. 




(a) (b) 



Figure 7: (a) Example of the cluster consisting of connected two polymers that 
share only one bond, (b) Example of the cluster consisting of disconnected two 
polymers that share only one site. 
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Figure 8: Example of the cluster consisting of connected four polymers that 
surround a plaquette and each neighboring pair of polymers share only one site. 
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Figure 9: L x x L y lattice (in this figure L x — 4 and L y = 5) with the fixed 
boundary condition for the ordered phase. The spins of the sites on the bound- 
ary (open circles) are fixed to a single value among 1, 2, • • • , q. 
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Figure 10: (a) Example of the configuration X of bonds (solid lines) for the 
free boundary condition, (b) The configuration X of removed bonds (crosses) 
for the fixed boundary condition with X <~ X. 
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q=b q=l y=8 5=10 




(2) 

Figure 12: The behavior of the ratio of the second energy cumulants F d to x. 



9 =20 




Figure 13: The behavior of the ratio of the absolute value of the third energy 



cumulants |i^ 3 } | to x 5 / 2 . 




Figure 14: The behavior of the ratio of the fourth energy cumulants to 
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Figure 15: The behavior of the ratio of the absolute value of the fifth energy 



cumulants \F^\ to x 11 / 2 . 
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Figure 16: The behavior of the ratio of the sixth energy cumulants F^l to x 7 . 



q=5 9=6 q=l ,7=8 g=10 



4=20 




0.2 0.4 0.6 0.8 1 1.2 1.4 

e 

Figure 17: The behavior of the comibined quantity of Eq.(21) in the disordered 
phase. 
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Figure 18: The behavior of the comibincd quantity of Eq.(21) in the ordered 
phase. 
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Tabic 1.1: The large-g expansion coefficients aJJ for the n-th energy cumulants 
of the free energy at /3 = (3 t in the disordered phase (n = 1-4). 



m 




(2) 


(3) 


(4) 



1 


2 


2 


2 


2 


2 


4 


14 


58 


242 


3 


-2 


26 


338 


3026 


4 


2 


118 


2474 


37186 


5 


-6 


250 


11490 


291946 


6 


4 


894 


55506 


2068914 


7 


-16 


1936 


219404 


12144136 


8 


6 


6160 


876384 


65992348 


9 


-38 


13538 


3127818 


322997186 


10 




39774 


11161914 


1494569130 


11 


-76 


88360 


37125056 


6470376688 


12 


-56 


245188 


122847224 


26851465312 


13 


-96 


547468 


387657840 


106378391164 


14 


-348 


1457976 


1214792028 


407805474096 


15 


156 


3264012 


3679570740 


1510368315660 


16 


-1634 


8410284 


11061186460 


5448643138308 


17 


1946 


18868858 


32420680250 


19136730128506 


18 


-6852 


47391870 


94306690802 


65779858359354 


19 


10744 


106180532 


269004235012 


221297704673276 


20 


-27004 


261607968 


761823806996 


731233034945892 


21 


48492 


586199668 


2124082332388 


2373735516063508 


22 


-102226 


1415498882 






23 


199386 


3174300542 
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Tabic 1.2: The large-g expansion coefficients aJJ for the n-th energy cumulants 
of the free energy at /3 = (3 t in the disordered phase (n = 5, 6). 



m 


5 

Mm 


(6) 
dm 



1 


2 


2 


2 


994 


4034 


3 


24338 


186626 


4 


488162 


5995138 


5 


5985042 


110362330 


6 


60977034 


1583071794 


7 


503300684 


17852715736 


8 


3668770296 


171621163420 


9 


23690762442 


1432076136578 


10 


140396232786 


10720828033434 


11 


768231161504 


72985931413600 


12 


3949989951344 


459566465021728 


13 


19184686377744 


2701568006995708 


14 


88891085678652 


14974796851686096 


15 


394525773133188 


78777820375010652 


16 


1687573606453276 


395842580467193364 


17 


6978567585546266 


1908650742049266298 


18 


28012258546890218 


8869906520648194890 


19 


109405587397512484 


39863282911219604012 


20 


416940286121010524 


173802434948834045508 


21 


1553287387425965092 


737013278292605504308 
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Table 2.1: The largc-g expansion coefficients bm for the n-th energy cumulants 
of the free energy at /3 = (3 t in the ordered phase (n = 1-4). 



m 


, i 

h y ' 


,2 
h y ' 


, 3 

h y ' 


h y ' 



1 


2 








2 


—4 


16 


—64 


256 


3 


2 


34 


-430 


3778 


4 


-2 


114 


-2654 


41778 


5 


6 


254 


-12186 


322670 


6 


-4 


882 


-57018 


2210982 


7 


16 


1944 


-224732 


12819264 


8 


-6 


6128 


-888024 


68657204 


9 


38 


13550 


-3164682 


333583598 


10 




39698 


-11243178 


1532324246 


11 


76 


88360 


-37363472 


6604807168 


12 


56 


245036 


-123377384 


27298396784 


13 


96 


547356 


-389128512 


107855738700 


14 


348 


1457784 


-1218076500 


412466192928 


15 


-156 


3263316 


-3688318020 


1524965526804 


16 


1634 


8410596 


-11080768444 


5492850472332 


17 


-1946 


18865590 


-32471142890 


19269581858838 


18 


6852 


47395762 


-94419894146 


66169209300214 


19 


-10744 


106166828 


-269288597908 


222430064188868 


20 


27004 


261629456 


-762460849076 


734462791941548 


21 


-48492 


586145660 


-2125652044660 


2382881224028156 


22 


102228 


1415594740 






23 


-199212 


3174081000 
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Table 2.2: The large-g expansion coefficients bm for the n-th energy cumulants 
of the free energy at /3 = (3 t in the ordered phase (n = 5, 6). 



m 


, (5) 


, 6 



1 

2 


-1024 


4096 


3 


-29518 


218914 


4 


-556382 


6813714 


5 


-6773994 


126069374 


6 


-67122114 


1774583142 


7 


-546094604 


19774354944 


8 


-3918393456 


187361651588 


9 


-25037212842 


1545876302510 


10 


-146961943266 


11451708807878 


11 


-798499755344 


77296110810160 


12 


-4080741048224 


483066658347536 


13 


-19726182681504 


2822025943834956 


14 


-91033421848212 


15558449192797824 


15 


-402728474968788 


81476650613568516 


16 


-1717926911825116 


407801271671111676 


17 


-7087982430811466 


1959732852092226390 


18 


-28396321463796218 


9080903533904186662 


19 


-110725077538793524 


40709407799692233428 


20 


-421377562607301884 


177104383529723970956 


21 


-1567944342623768692 


749588012269824549980 
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Table 3: Estimates of the specific heat Cd for the disordered phase in comparison 
with the low-temperature series and the Monte Carlo simulations. 



q 


large- q 


low-temp. 


Monte Carlo 


5 


2889.1(41) 






6 


205.930(55) 






7 


68.7370(52) 


68.01(47) 




8 


36.9335(11) 


30.17(6) 




9 


24.58756(35) 


20.87(5) 




10 


18.38542(17) 


16.294(34) 


18.437(40) 


15 


8.6540342(98) 


8.388(4) 


8.6507(57) 


20 


6.132159678(47) 


6.0645(69) 


6.1326(4) 


30 


4.29899341467(73) 


4.2919(51) 





Table 4: Estimates of the specific heat C for the ordered phase in comparison 
with the low-tcmpcraturc series and the Monte Carlo simulations. 



q 


large- q 


low-temp. 


Monte Carlo 


5 


2885.8(34) 






6 


205.780(32) 






7 


68.5128(22) 


52.98(61) 




8 


36.62347(31) 


28.47(47) 




9 


24.203436(65) 


20.03(27) 




10 


17.937801(18) 


15.579(83) 


17.989(40) 


15 


7.99645871(22) 


7.721(28) 


7.999(3) 


20 


5.360768767(13) 


5.2913(53) 


5.3612(4) 


30 


3.41289525542(28) 


3.4012(10) 
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Table 5: Estimates of the energy cumulants F\ n for the disordered phase (n = 
3-6). The values in the brackets are the data of the Monte Carlo simulations. 



q 


r d 


r d 


5 


2.1597(27) xlO y 


7.4380(99) xlO 15 


6 


2.04752(87) xlO 6 


1.07367(47) xlO 11 


7 


9.8404(16) xlO 4 


8.2877(14) xlO 8 


8 


16420.4(11) 


4.69139(37) xlO 7 


9 


4880.80(16) 


6.69421(26) xlO 6 


10 


2002.017(32) 


1.600281(33) xlO 6 




[2015(26)] 


[1.583(64) xlO 6 ] 


15 


176.78507(17) 


32323.501(59) 




[176.01(76)] 


[2.93(20) xlO 4 ] 


20 


54.8121419(72) 


4904.2443(16) 




[54.7(4)] 


[4905(89)] 


30 


15.88750157(25) 


667.244592(17) 



q 


P 5 
b d 


b d 


5 


4.90(39) xlO 22 


4.9(15) xlO 29 


6 


1.092(30) xlO 16 


1.68(38) xlO 21 


7 


1.363(14) xlO 13 


3.44(52) xlO 17 


8 


2.625(11) xlO 11 


2.30(20) xlO 15 


9 


1.8024(36) xlO 10 


7.67(37) xlO 13 


10 


2.5155(25) xlO 9 


6.29(16) xlO 12 




[2.40(26) xlO 9 ] 


[5.8(18) xlO 12 ] 


15 


1.168612(78) xlO 7 


6.779(13) xlO 9 




[1.170(59) xlO 7 ] 


[7.2(11) xlO 9 ] 


20 


870078.5(88) 


2.47889(69) xlO 8 




[8.42(32) xlO 5 ] 


[2.33(16) xlO 8 ] 


30 


55619.490(41) 


7.44878(13) xlO 6 
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Table 6: Absolute value of the estimates of the energy cumulants Fo for the 
ordered phase (n = 3-6). The values in the brackets are the data of the Monte 
Carlo simulations. 



q 




F^ 


5 


2.1577(84) xlO 9 


7.368(33) xl0 ib 


6 


2.0467(28) xlO 6 


1.0698(22) xlO 11 


7 


9.8469(57) xlO 4 


8.2777(91) xlO 8 


8 


16460.8(46) 


4.6939(32) xlO 7 


9 


4905.23(73) 


6.7098(35) xlO 6 


10 


2018.31(17) 


1.6071(11) xlO 6 




[2031(26)] 


[1.55(22) xlO 6 ] 


15 


181.4445(19) 


32878.80(86) 
[2.93(20) xlO 4 ] 




[180.67(76)] 


20 


57.05595(23) 


5054.1460(86) 




[57.09(29)] 


[4.79(22) xlO 3 ] 


30 


16.7446548(22) 


702.788233(56) 



q 


F^ 




5 


5.07(21) xlO 22 


5.15(46) xlO 29 


6 


1.105(13) xlO 16 


1.775(60) xlO 21 


7 


1.3691(65) xlO 13 


3.586(52) xlO 17 


8 


2.6335(57) xlO 11 


2.356(16) xlO 15 


9 


1.8084(19) xlO 10 


7.799(27) xlO 13 


10 


2.5268(14) xlO 9 


6.367(12) xlO 12 




[1.98(51) xlO 9 ] 


[2.8(14) xlO 12 ] 


15 


1.185486(57) xlO 7 


6.8744(12) xlO 9 




[0.92(15) xlO 7 ] 


[4.5(13) xlO 9 ] 


20 


892352.53(68) 


2.538205(72) xlO 8 




[7.37(70) xlO 5 ] 


[1.38(26) xlO 8 ] 


30 


58121.979(30) 


7.760343(17) xlO 6 
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